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Abstract 

In this paper we investigate the statistical properties of gravitational lenses 
for models in which a cosmological term decreases with time as A oc a~"^, where 
a is the scale factor and m is a parameter (0 < m < 3). We show that for given 
low values of the present matter density parameter > there is a wide range 
of values for m for which the lensing rate is significantly smaller than that in 
cosmological constant (A) models. We also show that models with low i^mo 
and m ~ 2 have high likelihood to reproduce the observed lens statistics in the 
HST snapshot survey. 
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1 Introduction 



In the last five years or so, the statistics of gravitational lensing (Turner, Ostriker 
and Gott III 1984, hereafter TOG) has proven to be a powerful tool in constraining 
models of the universe, especially those with a cosmological constant (A). Cosmologies 
with a A-term have a long history and are now receiving considerable attention (see 
Carroll, Press and Turner 1992 for review). Spatially flat cosmological models with a 
cosmological constant have been suggested (Peebles 1984; Turner, Steigman & Krauss 
1984) as a way to reconcile inflation with dynamical analyses on scales of ~ lOhr^ 
Mpc, that indicate a value for the density parameter flo ~ 0.1 to 0.3 (Peebles 1993). 
A cosmological constant also alters the transfer function for the density perturbations 
giving more power in the perturbation spectrum at large scales (as compared with 
standard CDM) in accordance with observations (Efstathiou et al. 1990; Lahav et 
al. 1991; Kofman et al. 1993). Besides, if the present value of the Hubble parameter 
is high, as indicated by some recent observations (Pierce et. al. 1994; Freedman et. 
al. 1994), a cosmological term will be the only way to get a theoretical age for a 
flat universe in accordance with current age estimates for globular clusters (Chaboyer 
1994). 

The idea that light could be focused by the gravitational lens effect was flrst 
suggested by Lodge (1919) near the beginning of the century. For several decades 
the subject of gravitational lensing had a quite slow development, but recently it 
started to become one of the most active research area in astrophysics and cosmology. 
There are several reasons for the current interest in gravitational lensing. On the 
cosmological side, after the works of Rcfsdal (1964) and Press & Gunn (1973), it was 
realized that cosmological parameters could be probed by the gravitational lensing 
effect. In the beginning it seemed that lensing properties were too insensitive and 
would only distinguish extreme cosmological models. Later on. Turner (1990) and 
Fukugita, Futamase & Kasai (1990) showed that a non-zero cosmological constant 
could signiflcantly affect the statistics of gravitational lenses, especially in a low- 
density universe. 

However, there are uncertainties in the study of the statistics of gravitational 
lensing (Mao 1991; Fukugita, Futamase, Kasai & Turner 1992, hereafter FFKT). 
For example, the lens effect depends considerably on how the mass is distributed 
in the lensing galaxy. Hinshaw and Krauss (1987) showed that the introduction of 
a core in the isothermal sphere galaxy model (non-singular isothermal sphere) can 
signiflcantly modify the statistical lensing properties. Another issue is what distance 
formula should be used: the angular diameter distance or the Dyer-Roeder distance? 
Related to this question is the kind of statistics to be applied (Ehlers & Schneider 
1986; FFKT). We should also mention the important effect of magniflcation bias and 
other selection effects such as angular resolution, galaxy evolution and merging on 
lensing probabilities (TOG; Fukugita & Turner 1991, hereafter FT; Mao 1991; Mao 
& Kochaneck 1993; Rix, Maoz, Turner & Fukugita 1994). 
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In spite of the uncertainties, the calculated rate of lensing in A-flat models, when 
confronted with the existing lensing observations, indicates that models with density 
parameter (Qmo) close to unity are most likely. For instance, Maoz & Rix (1993) 
claim that at present we should have Q\ ~ 0.7 . So, it is becoming more and more 
difficult to make the dynamical estimates for Q on scales ~ 10h~^ Mpc compatible 
with a fiat cosmological model with A 7^ 0. It should be pointed out, however, that 
the lensing frequency of quasar images is considerably reduced if early-type galaxies 
{z ~ 0.5) were dusty (Fukugita & Peebles (1994)). 

The results above seem to favor open FRW models. There are however other 
possibilities. For example, Ratra and Quillcn (1992) showed that, for a wide range of 
parameters, the predicted lensing rate is considerably reduced in a class of fiat models 
(Peebles & Ratra 1988, Ratra & Peebles 1988) in which a scalar field plays the role 
of an effective cosmological "constant" that decreases with time. Some other models 
with a decreasing cosmological term were also proposed ((Ozcr & Taha 1987a,b; 
Freese, Adams, Frieman & Mottola 1987; Chen & Wu 1990; Abdel-Rahman 1992; 
Carvalho, Lima & Waga 1992; Silveira & Waga 1994) and it would be interesting 
to know if they also predict a lower lensing rate. We can argue that we should 
expect a positive answer to this question. The reason is that usually in a varying 
A cosmological model, the distance to an object with redshift z is smaller than the 
distance to the same object in a constant A model with the same Qmo- So, the 
probability that light coming from the object is affected by a foreground galaxy is 
reduced in a decaying A cosmology. However this is only a qualitative argument, 
and it is clear that a quantitative treatment is necessary if we want to put limits on 
parameters of the models. 

In this paper we address the above question to the special class of models proposed 
by Silveira and Waga (1994) in which a cosmological term decreases with time as 
A cx a"'", where a is the scale factor and < m < 3 is a constant. We show that 
these models also admit a large set of parameters for which the predicted lensing rate 
is much lower than that obtained in a constant A model with the same low value 
of flmo- The paper is organized as follows: In section (2) the assumptions and basic 
equations of our models are presented. We exhibit expressions for two sorts of distance 
that we shall use, and discuss the corresponding statistics associated with them. In 
section (3) we model galaxies by the isothermal sphere density profile and obtain the 
predicted lensing probabilities and the distribution of image angular separation for 
some typical models previously chosen. In section (4) we compare the predictions of 
the models to observations and stress our main conclusions in section (5). 

2 Decaying vacuum cosmological models — distance 
and optical depth formulas 

Following Silveira and Waga we assume that the cosmic fluid is a non-interacting 



2 



mixture of non-relativistic matter and radiation. The cosmological term is assumed 
to be a time dependent quantity, 

A = spa-"", (1) 

where /9 > is a constant and the factor 3 was inserted for mathematical conve- 
nience. We also assume that the vacuum decays only into relativistic particles such 
that matter is conserved {pm oc a^^ ). As shown by Silveira and Waga (1994), the 
radiation energy density has two parts; one conserved {QroHo^{(io/ay) and a second 
one, ( g^J(^^^^) Q~"')i that appears due to particle creation by the decaying vacuum. 
Here ao is the present value of the scale factor, Hq — 100 h km Mpc^^ is the 
present value of the Hubble parameter {h ~ 0.5 — 1) and Qro = 4.3 x 10^^/i~^ stands 
for the present value of the conserved radiation density parameter. In the following, 
subscripts will always indicate present values. 

The Einstein equations for the models we are considering reduce to 

[l] =^-o^o(^) +^.oH',(^^) (2) 

and 

- - -^^moH, + -^—n^oH, (^-J , (3) 

where flrno is the matter density parameter, Q^o = — ^^=0 = 

Since we are mainly interested in the lensing properties of the models, only recent 
epochs have to be considered (2; ~ 5). This justifies neglecting the conserved radiation 
energy density on the right hand side of (2) and (3). To have some grounds of 
comparison we have included the curvature term in (2) and will also consider the 
open FRW (A; = — 1) model. 

The equations above are quite general and apply for a broad spectrum of models. 
Let us first consider the A; = case. For instance, if m = the usual flat FRW model 
with a cosmological constant is recovered, whilst if m = 2 the above equations (with 
^fco = 0) are formally the same as those of the open FRW model. The same equations 
also appear in some string dominated cosmologies (Vilenkin 1984). Further, we would 
get the same equations if we had considered, besides matter , an exotic x-fluid with 
equation of state, Px — — ^)Px- Cosmologies having a fluid with this behavior were 
investigated by Fry (1985), Sahni, Feldman and Stebbins (1992), Feldman and Evrard 
(1993) and more recently by Martel (1995). We remark that all we shall discuss here 
also applies for these models. We should also mention that in the limit Pm » P4>, 
the scalar field model analyzed by Ratra and Quillen has the same behavior as the 
one proposed here. This can be seen easily if we relate their parameter a with m 
as, a — It is clear however, that the models are different (unless m = 0) if 

Qrno ~ ^xo or duriug the x-component (vacTium) dominated era when, in fact, all the 
lensing properties we shall discuss are important. 
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In this paper we shall compare the following models: 
Case A: m — 0, Q^o — 1 and A; = (Einstein-de Sitter model). 
Case B: m = 2, Q^q — 0.2 and k — 0. 

Case C: m = 0, flmo = 0.2 and k = (Friedman-Lcmaitre model). 
Case D: f^^o = 0.2, fi^o = and A; = -1 (Open FRW model). 

For the sake of completeness we have included case D in our analysis. It will be 
interesting to compare it with case B, which has similar field equations but has flat 
spatial sections. We also analyzed the case m — 1 and it turned out that it always 
has behavior between cases B and C, so we decided not to explicitly include it in our 

discussion. 

For the flat models A, B and C, the angular diameter distance, ds{zL, zs) , between 
two objects, one with redshift and the other with zs is given by, 

ds{zL, zs) = — ^ / , =• 4 

^ + ^sJz, ^O^o(l + |/)3+(l-f^mo)(l + ?/)"* 

Equation (4) can be expressed in terms of the hypergeometric functions F{a, b; c, z) 
as 

[1 + Zs)V^l-mO 

( ^ ^ 7 — 2m 1 — ^mo \m-3\ _ 

^VT+^ 4' 6 -2m' 6 -2m' Qmo ^ ^ ''^ ^ 

^yl + Zs 2 D — 2m D — 2m y^o 

For some special values of m, the hypergeometric function in (5) can be reduced to 

elementary functions. This can be done, for instance, for m = 2 if we use the relation, 
F(l/2, 1/2; 3/2, —x'^) = sinh~^(x). For the Einstein-de Sitter model equation (4) 
can easily be integrated giving 

dsizL, zs) = [(1 + z,)-'/' - (1 + zs)-'/'] . (6) 

1 + Zs ^ 

In fact, we can obtain (6) from (5), by observing that in the limit Jl^o 1 the 
hypergeometric function also goes to unity. 

In the case of open models, the angular diameter distance can be expressed as 
(FFKT) 

2cH~' 

ds{zL, Zs) = „2 , TT:— r^{(2 - ^2„^o + ^mOZs) X 

^ mo(l + ^i:) (l + zs)^ 

\Jl-\- f^mO^L - (2 - f^mO + ^mO^L)\/l + ^mO^s}- (7) 
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The differential probability, dr, that a line of sight intersects a galaxy at redshift 
zl in the interval dzi from a population with number density rii is (TOG; Peebles 
1993) 

dr — —niac- — dzL, (8) 

dZL 

where from (2) we have 

I = = (^^..0 (1 + zf + (1 + zr - ^ko (1 + zf)-'^' . (9) 

The cross section (cr) in (8) is given by 

a^nal, (10) 

where the effective radius of the lens, that is^ df^f IS the maximum distance of 

the lens from the optical axes for which multiple image is possible. 

The total optical depth (r) is obtained by integrating dr along the line of sight 
from to zs, that is 

fZS fZS dt 

T — dr — — niGC- — dzL- (11) 
Jo Jo dzL 

In the angular diameter distance definition it is assumed that the matter in the 
universe is homogeneously distributed. However the gravitational lens effect will not 
occur in a smooth imiverse. Only if matter is clumped, as in the real universe, can 
this effect take place. A distance formula that takes matter clumping into account 
was proposed by Dyer and Roeder (1972, 1973) and is known as the Dyer-Roeder 
(DR) distance. Here we will consider two extreme cases. We have already discussed 
the first one in which the smoothness parameter (a), where < a < 1, is equal to 
one (filled beam DR distance or angular diameter distance). The other extreme case, 
Q! = 0, is called the DR empty beam distance (Schneider, Ehlers and Falco 1992). 

For the models under consideration, the empty beam DR distance is given by ( 
FFKT) 



dDR{zL,zs) = cHq {1 + zl)x 



y) dy 



=(12) 



+ yf + (1 - - J^xo)(i + y? + 5^x0(1 + yy 

Notice that for the same ilmo, flat models (Qmo + = 1) with m — 2 and open 
models {VL^q = 0) have the same empty beam distance. For the open ( just make m — 2 
in (13) ) and flat models, equation (12) can be expressed in terms of hypergeometric 
functions as 
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, , , 2cH^\1 + zl) 

aDR[ZL,Zs) = ^ — X 

Again, in the limit 1 the hypergeometric function goes to unity and (13) 

simphfies to 

dnnizL, zs) = ^^^^^i^((l + z,)-'l' - (1 + zsY'''). (14) 

In obtaining the probabihty of multiple images in (11), we considered a random 
line of sight to the source at zs, calculated the expected number of lenses {dr) in 
the redshift interval dz^ around z^, and then integrated dr from to zs- Ehlcrs 
and Schneider (1986) observed that in a self-consistent treatment of probabilities in a 
clumpy universe, the random variable should be the position of the source on a sphere 
at Zs (and not the line of sight to the source). They then proposed a new derivation for 
the optical depth that is called the ES probabihty. The Ehlers-Schneider differential 
probability (drEs) can be expressed as (FFKT) 

where in dr (given by equation (8)), the empty beam distance should be used. By 
integrating (15) from to zs we obtain the total ES optical depth. 

In Figures la and lb we present the quantity Dqs/cHq^ for the filled and empty 
beam distances. We also show in Figures 2a and 26 , also for both distances, 
the combination DqlDls/ {DqscHq^) that appears, through acr, in the differential 
probability formulas. We are following the TOG and FFKT notation, such that 
Dls = d{zL, Zs), Dos = d{0, zs) and Dql = d{0, zl). 



3 The isothermal sphere galaxy model 

Let us now consider the isothermal sphere model for galaxies. This model is 

characterized by two parameters, namely, the core radius (rc) and the one-component 
velocity dispersion (o"||). Following Hinshaw and Krauss (1987), we assume the lens 
galaxy to have the following density profile 
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The surface mass density of the lens on the lens plane is given by (Bourassa & 
Kantowski 1975) 



o 

E(a)=2/ P(r)^^= I ^ , (17) 



Vr2 - a2 2Gx/^^T77^ 
and the projected mass interior to the impact parameter b is, 



/•o Tran r / 

M(6) = 27r / aE(a)cia = — ^ a/^^ + - 
^0 G L 

The bending angle is 



(18) 



« = ^^ = «o ^ , (19) 

where ao — Att (^^^ 1.8" (cr||/250 kms'^Y denotes the constant bending angle for 
the singular isothermal sphere case (SIS), obtained by taking the limit Tc — > 0. 
By using simple geometry it is easy to see from Figure 3 that 

l + b^ — a, (20) 

^os 

where I is the distance from the lensing galaxy to the unperturbed line of sight. It 
follows from (19) and (20) that, if Tc = (SIS case), the maximum value of I for 
multiple images (ocr) is given by 

«cr(0) = — ao- (21) 

J-'os 

If Tc 7^ (NSIS), by substituting (19) in (20) we get the cubic equation, 

b^ + 2lb^ + {l^ + 2reae,(0) - al{0))b + 2lrca^{0) = 0. (22) 

The number of real and distinct solutions of (22) depends on the sign of its dis- 
criminant. Hinshaw and Krauss (1987) showed that in order to produce multiple 
images, the lens maximum distance from the unperturbed line of sight should have 
the following expression 



(23) 



where (3 — rc/acr{0). Further, from the multiple image diagram (Young et. al. 
1980, Blandford and Kochanek 1987), they also showed that multiple images are only 
possible if /3 < 1/2. So, the cross section for NSIS is cr = for [3 > 1/2 and a = vra^^ 
(with Qcr given by equation (23)) if < 1/2. In fact, if /? < 1/2, instead of two 
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(as in SIS), three images are predicted, in agreement with the odd number of images 
theorem vahd for symmetric (non-singular) lenses (Subramanian & Cowling 1986). 
It can also be shown (Hinshaw and Krauss; see also Hinshaw 1988 for more details) 
that the angular image separation between the two outer images is given by 



e 



2a, 



Do 



L 



(24) 



1/2 



Notice that if /? = the usual result 9 = Stt^^ (^)' , valid for SIS lenses, is recov- 
ered. 

In the SIS case, by using equation (21) and assuming conserved comoving number 
density of lenses ( rii = no{l + z)^ ), equation (11), can be analytically integrated. By 
using standard distance and statistics we obtain for flat models. 



Here, 



T{zs)^^{ds{0,zs){l + zs)f. (25) 

/ = ^ < noaf^ > (26) 

measures the effectiveness of the lens in producing multiple images (TOG). An ana- 
lytic expression for r in the case of filled beam distance and standard statistics can 
also be obtained for k ^ (see Gott, Park and Lec 1989). 

Following FT, we consider the existence of 3 species of galaxies (E, SO and S) and 
assume a Schechter form for the luminosity function, 

= 0*(^)"exp(-L/L*)^, (27) 

where 0* = (1.56 ± 0.4) x IQ-'^h^Mpc-^ (Efstathiou, Ellis & Peterson, 1988) is a 
galaxy number density and a = —1.1 ± 0.1 (see FT). The morphological composition 
is assumed to be £■ : SO : 5" = 12 : 19 : 69. We assume, in addition, the relationship, 
{^) — {■^)'^, between galaxy luminosity and velocity dispersion. Here the exponent 

7 is, 7 = 4 for E/SO ( Fabcr and Jackson 1976) and 7 = 2.6 for S galaxies (TuUy and 
Fisher 1977). Substitution of (27) in (26) leads to, 

1 flTT^ 4 

c-Hq 7 

where r(x) is the Gamma function and a*^ is the velocity dispersion corresponding to 
the characteristic luminosity L*. 

Fukugita and Turner estimated crjl to be: 

2251:11 /or E, 
(7* = { 206tll for SO, (29) 
1441^3 for S. 
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To take into account dark massive halos, they follow TOG and also adopted a (3/2)^/^ 
correction factor for the velocity dispersions for E/SO galaxies. With the above 
numbers we find /e = 0.018 ± 0.009, fso = 0.020 ± 0.011, and fs = 0.007 ± 0.003 
(total effectiveness parameter, / = 0.045 ± 0.023). More recently Kochanek (1993) 
argued that the ratio of dark matter dispersion velocity to that of luminous matter 
should be in the range 0.9 — 1.05 and suggested that the (3/2)^/^ correction factor 
should not be considered. Without this the factors /e/so a-re 2.25 smaller, that is, 
fE = 0.008 ± 0.004, fso = 0.009 ± 0.005 and we get / = 0.024 ± 0.012. We shall 
consider in the next section these two cases when comparing the predictions of the 
models with observations. 

In Figures 4a and 46 the normalized optical depth (r//) , for the four models in the 
SIS case, is displayed for filled beam distance and standard statistics (Figure 4a) and 
for empty beam distance and Ehlers-Schneider statistics (Figure 46). We also obtained 
the optical depth for the NSIS case. As discussed before, in this case we should use 
the appropriate cross section with a^ given by (23). For an analytic expression of 
the NSIS optical depth in the standard case see Krauss and White (1992). In Figure 
5a is displayed the NSIS normahzed optical depth for filled beam distance (standard 
statistics). In Figure 56 the same quantity is displayed for the empty beam distance 
(Ehlers-Schneider statistics). A constant value for the core radius, Tc = 0.5h''^kpc 
and a velocity dispersion a*^ = 144 km/ s are assumed for all the models. 

The present available data do not allow very good estimates for the core radius of 
galaxies. In fact, Tc seems to vary a lot even among galaxies with the same morphol- 
ogy. However there is some evidence that most spirals have large core radius {tc ~ 0.5 
/cpc), although there are also indications that ~ 10% of them have very small cores 
(FFKT and references therein). So, in view of the lack of more precise information 
we assume that 90% of spiral galaxies have a constant core Tc = 0.5h^^kpc and that 
the remaining 10% are well described by SIS. In fact these assumptions are enough to 
practically reduce the contribution of S galaxies to the optical depth to only 10% of its 
SIS value. Actually in case C (cosmological constant), remains another ~ 1% effect. 
We can understand this small difference by observing that the quantity DolDls/ Dqs 
for fixed redshift is higher in case C (see Figure 2). This means that the parameter 
13 tends to be smaller in case C and explains why the effect should be less important 
in this case. In any case, we confirm the conclusion obtained by Krauss and White, 
that spiral galaxies have a very small effect on lensing frequencies even if we do not 
include the ^3/2 factor in the E/SO dispersion velocities. 

While spirals usually have large core radius, E and SO galaxies arc believed to 
have smaller ones. Most of the analyses of elliptical galaxies are based on slightly 
different relations between core radius and velocity dispersion (luminosity) that are 
usually derived by using Lauer's (1985) study of nearby elliptical galaxies. Laucr 
found that 14 galaxies (from a total number of 42) have resolved cores (SOpc ~ Tc ~ 
400pc), 23 were unresolved and 5 had marginally resolved cores. By fitting Lauer's 
data for the resolved cores, Krauss and White obtained a relation between E core 
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radius and luminosity. By assuming that relation to be valid for all E galaxies they 
obtained a suppression factor ~ 0.4 for the Einstein-de Sitter model and ~ 0.63 for 
the cosmological constant dominated universe. Fukugita et al. (FFKT) assumed 
that 1/3 of E galaxies are well described by the relation they obtained from Lauer's 
study. They also assumed that another 1/3 have core radius given by multiplying 
that relation by 1/3 and that the remaining 1/3 have Vc ~ lOpc. With this model 
they obtained a suppression factor equal to 0.65. They also claim having changed 
their assumptions in a reasonable way and always getting numbers between 0.5 and 
0.7. We analysed this effect more quantitatively and reached similar results. 

We also analysed the core radius effect in the case of empty beam and ES-statistics. 
We found that the suppression is higher in this case. The reason is that acr(O) for 
the empty beam case is smaller and this implies that (3 is higher, thus increasing 
the suppression effect. Actually, for spiral galaxies the effect can be very high. For 
instance, from Figures 4a and 46 it is clear that Tsisi"^, filled beam)/rs'/5(z, empty 
beam) is less than one. This means that we should expect a higher frequency of 
lensed quasars in the open beam case. However by looking at Figures 5a and 56 
we immediately realize that for the special choice of the parameters T^sisi'^', filled 
beam)/T7V5/5(z, empty beam) is higher than one and the opposite would be expected. 
In fact this occurred because we considered in our example a typical spiral galaxy 
with relatively high core radius and small velocity dispersion. For E/SO galaxies 
we should expect this effect not to be so conspicuous and, in fact, under reasonable 
assumptions we obtained a suppression factor that is only ~ 10% higher than that 
in the filled beam case. In the next section, when comparing model predictions to 
observations, we will take this small difference into consideration. In order to simplify 
calculations, we will follow FFKT and assume, for filled beam distance, a constant 
core effect suppression factor equal to 0.65 for E/SO galaxies. In the case of the 
empty beam we will consider a suppression factor equal to 0.60. 

By using the filled beam distance, standard statistics and the SIS profile it can be 
shown (see FT and FFKT) that the normafized image angular separation distribution 
for a source at zs is given by 



dP _ 1 r^s 



7/2 



r(a + i + 



/ 



T{zs} Jo 



OS 



D 



LS 



\ 



8tt 



i(«+2) 




exp 



(30) 



In Figure 6a the predicted image separation distribution for a source at redshift 
Zs — 2.2 (the quasar average redshift in the HST snapshot survey) is shown in two 
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cases for E/SO galaxies. We took a = —1.1, and to simplify the computation we 
considered the same velocity for E and SO galaxies, that is we chose an average 
velocity, a*^ — 213 x km/s in case (i) and a*^ — 213 km/s in case (ii). The same 
quantity for both cases is also displayed in Figure 6b for the empty beam distance 
and Ehlers-Schneider statistics. In this case we took into account equation (15) in the 
definition of the differential optical depth. Figure 6a shows that for fiat models and 
filled beam distance the image separation distribution is independent of m and flmo- 
In fact, as remarked by FFKT, it is also independent of zs- For empty beam distance 
the degeneracy of fiat models is broken and we can observe a shift of the distribution 
to smaller values of angular separation. It is also clear that in both cases increasing 
the velocity dispersion increases the probability of larger image separation. 

4 Comparison with observations 

In this section we shall compare the predicted number of multiple images (and their 
angular image separation) for the models presented in section 2 with the observational 
results of the Hubble Space Telescope Snapshot Survey (Maoz et al. 1993). In the last 
report of the survey, Maoz et al. announced the existence of six lens candidates from 
a sample of 502 quasars. In fact, two of the lensed candidates have unexpectedly large 
image separation. There is evidence that one of them was produced not by a single 
galaxy but by a more complex system, and it is not clear that the other one is really a 
multiple image of one single object (Maoz and Rix and references therein) . Following 
the current interpretation, we shall not include these two cases in our analysis. The 
four remaining candidates have the following redshift and image separation {z,9): 
(3.8, 0.47"), (2.55, 1.22"), (1.72, 2.0") and (2.72, 2.2"). 

The expected number of lensed quasars in the survey {Ne), for each model, is 
obtained by computing the quantity 

Nq 

NE = Y.r{zi) (31) 

i=l 

where the sum is over the Nq quasars rcdshifts of the survey. In calculating the 
expected number of multiple images, besides the core radius effect that we discussed 
in the last section, we also took into account two other corrections to r, namely, the 
angular selection effect ( xO.95 for E and SO) and magnification bias (x9.1) (see 
Fukugita and Peebles). 

In Table 1 we display, for each model, the predicted number of lensed quasars for 
the HST snapshot survey. The numbers were obtained for filled beam and standard 
statistics as well as for empty beam and Ehlers-Schneider statistics. We show results 
for the cases in which we have and have not included the (3/2)^/^ factor in velocity 
dispersion. By assuming Poisson statistics we also display, for all cases, the probability 
of detecting the observed number (four) of lensed quasars. Notice that in case C the 
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final results depend considerably on the assumptions performed. For instance, if we 
take into account the (3/2)^/^ factor, the predicted number of lensed quasars is too 
high in model C. However if we do not consider the correction the model cannot be 
ruled out, at least based only on the total number of images. From Table 1 we also 
see that model D is the least sensitive to the assumptions and seems to be in good 
agreement with observations. It is also clear that the lensing rate is a factor ~ 2 
smaller in model B than in model C and this corroborates the idea that usually decay 
A models have a lower lensing rate than the constant A ones. We should remark that 
this is not always true and some decaying A models can, in fact, predict very large 
lensing rates. 

Recently Kochanek (1993) pointed out that analyses based only on Poisson statis- 
tics give no weight to the different ways that a fixed number of multiple images is 
produced. He then suggested a maximum likelihood method that takes this into 
account. His technique is based on the following likelihood function, 

L=\{{l-pi)\{p,\{p,k. (32) 

i=l j=l k=l 

Here Nu is the number of unlensed quasars, A^^ is the number of lensed ones, Pi <^ 1 
is the probability that quasar i is lensed and p^fc is the configuration probability, that 
we shall consider as the probability that quasar k is lensed with the observed image 
separation. 

We applied Kochanek technique to the flat models. By expressing L as a function 
of the parameters m and Qmo we obtained the maximum of the likelihood function 
{Lmax) and formed the ratio / = L/ Lmax- It can be shown that with two parameters, 
the distribution of —2 In/ tends to a distribution with two degrees of freedom 
(Kendall & Stuart 1977, Eadie et al. 1971, Kochanek 1993). 

Contours of constant likelihood are plotted in Figure 7. Regions with larger like- 
lihood are represented by lighter shades. The maximum of the likelihood is indicated 
by a cross (-I-). We should remark that the figure is very broad and is displayed 
to give a qualitative view of the likelihoods. For the figure we considered the value 
of the velocity dispersion equal to 213 km/s. We observed that the results did not 
change appreciably when we increased the velocity dispersion to a*^ — 261 km/s. 
Some qualitative aspects can be inferred from the figure. For instance, if m is low 
(m ~ 0.5) regions with a lower value of A, let say, those with Vtmo ~ 0.4, have higher 
likelihood. However if m ~ 2 that constraint does not exist and models with low 
Q^o are in fact more likely. We observe that in the two-dimensional space of the 
parameters of Figure 7, the Einstein-de Sitter models is represented by two points. 
The first one is m = and VL^o = 1, while the second one is obtained by taking the 
limit m — * 3 and Vtmo — * 0. Figure 8 gives more quantitative information. We plot in 
it the 50%, 68%(1(t) and 95.4%(2(t) confidence levels for the likelihood ratio for the 
two parameters. 
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It is important to emphasize that our main goal in this paper was not to obtain 
constraints on possible values of a constant A. Our target in this work was rather 
to show that if A is not constant, the lensing constraints on fl\ are much weaker. 
In fact we found that regions with higher values of m and lower Qmo present larger 
likelihood. However, if we want to make some comparison with previous results, as 
for instance those of Kochanek (1993), we should fix the parameter m to the value 
m — (cosmological constant). In this case the likehhood peaks at Jl^o — 0.62 
but now we have that —2\nl is distributed like a distribution with one degree 
of freedom. In this case the constraints on A are stronger and, for instance, we get 
Qa ~ 0.8 at 90% confidence level. We remark that we have not obtained one of 
Kochanek's constraints, fl\ ~ 0.45, for the following reasons. First, as we discussed, 
we were only considering lensing by galaxies, so we have not included in our analysis 
the lens 0957+561. Second, we have used the SIS profile but we took into account 
a core radius suppression effect. Third, in our likelihood analysis we neglected the 
contribution of spiral galaxies. 

5 Discussion and conclusion 

As is well known, a sufficiently long period of inflation in the early universe gives a 
natural solution to the isotropy, flatness, and monopole problems. A nearly invariant 
primordial spectrum is also generated in most of the inflationary models of the uni- 
verse, a feature that seems to be in good agreement with observations. Inflationary 
models usually predict fltotai = 1- Nevertheless, observations on the scales 10 — 30 
Mpc, based on dynamical methods, indicate Q^o = 0.2 ± 0.1. Further, in flat models 
with Qjno — 1 only ii h < 0.59 is it possible to get theoretical ages in agreement 
with the lowest age estimates (to = llGyr) for globular clusters. However, most of 
the recent observations indicate higher values for the Hubble parameter. Besides, the 
standard CDM model {flmo = 1 and h = 0.5) when normalized to COBE predicts 
more power on small scales than is observed, and some of its variants such as HCDM 
seem to work well only if /i ~ 0.5 (Primack 1995). By assuming a non-zero cosmologi- 
cal constant all these problems can be solved at once while keeping the attractiveness 
of inflation. 

Cosmologists in general show an enormous resistance in accepting the idea of 
a non-zero cosmological constant for several reasons. First, because it is another 
parameter in the theory and from an aesthetic point of view this makes A-models 
less compelling. The second reason is that in order to dominate the dynamics of the 
universe only recently, this parameter should have a very small value (A ~ lO^^^cm^^) 
that is in fact 50 to 120 orders of magnitude below the estimate given by quantum 
fleld theory. If we assume a decaying cosmological term, this second problem is 
alleviated but we have to pay the price of introducing another parameter. In this 
paper we have also discussed a weakness on the observational side of constant-A 
models. It is the prediction of too high frequency of lensed quasars by models with 
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a large cosmo logical constant. Ratra and Quillen showed that the predicted lensing 
rate is considerably reduced in certain models in which a scalar field plays the role 
of an effective decaying A-term. In this paper we reached the same conclusion for 
models in which the cosmological term decreases with time as A oc a~™. We went one 
step further and also showed that for these models, lower values of Qmo and larger 
m have higher likelihood. Finally, we should mention that by increasing the value 
of the parameter m, and maintaining constant the value of flmo, the theoretical age 
predicted by the models we have considered decreases. So, it would be interesting to 
extend the likelihood analysis by also taking into account the age constraints. Further 
investigation in this direction is being carried out. 
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Table 1: Predicted number of lensed quasars for the HSTSS and model 
probabilities 

model A model B model C model D 

A; = A; = 0,m = 2 A; = 0,m = k ^ -1 

^mO — 1 ^mO — 0.2 ^mO — 0.2 ^mO — 0.2 



filled beam (with ^3/2) 


2.6 


4.7 


10.7 


4.1 


probability 


18% 


17% 


0.3% 


19% 


empty beam (with ^3/2) 


4.2 


7.1 


14 


4.6 


probability 


19% 


6% 


10-2% 


18% 


filled beam (without ^3/2) 


1.2 


2.2 


5.0 


1.9 


probability 


9% 


16% 


16% 


14% 


empty beam (without ^^3/2) 


1.9 


3.3 


6.8 


2.1 


probability 


14% 


20% 


7% 


15% 
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Figure Captions 



Figure 1: Dyer-Roeder distance as a function of the redshift; (a) filled beam and (b) 
empty beam. 

Figure 2: The combination DqlDls/ {DqscHq^) as a function of zl ; (a) filled beam 
distance and (b) empty beam distance. 

Figure 3: Lensing geometry. 

Figure 4: The normalized optical depth (r//) as a function of the redshift for the 
SIS case; (a) is for filled beam distance and standard statistics and (b) for empty 
beam distance and Ehlers-Schneider statistics. 

Figure 5: The normalized optical depth (r//) as a function of the redshift for 
the NSIS case; (a) is for filled beam distance and standard statistics and (b) for 
empty beam distance and Ehlers-Schneider statistics. For all the plots we assumed 
Vc — 0.5h~^kpc and a*^ — 144 km/s. 

Figure 6: Distribution of the angular image separation for a source at 2:5' = 2.2 in 
SIS case. We took a = —1.1, j — 4 and considered in case (i) a*^ = 261 km/s and in 
(ii) (Tj| = 213 km/s ; (a) filled beam and standard statistics and (b) empty beam and 
Ehlers-Schneider statistics. 

Figure 7: Contours of constant likelihood for fiat models are plotted in the flmo and 
m parameter space. Regions with larger likelihood are represented by lighter shades. 
The peak of the likehhood is indicated by the cross. 

Figure 8: Likelihood contours at 50%, 68%{la) and 95.4%(2(t) confidence levels for 
the two dimensional likelihood distribution I. 
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